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ABSTRACT 

Radiative efficiency, rj = L/Mc 2 , of a slowly accreting black hole is estimated 
using a two-temperature model of accretion. The radiative efficiency depends on 
the magnetic field strength near the Schwarzschild radius. For weak magnetic 
fields, /3~ L = B 2 /8iTp <^ 10" 3 , the low efficiency rj ~ 10~ 4 assumed in some 
theoretical models is achieved. For > 10~ 3 , a significant fraction of viscous 
heat is dissipated by electrons and radiated away resulting in 77 > 10~ 4 . At 
equipartition magnetic fields, fi^ 1 ~ 1, we estimate 77 ~ 10 _1 . 
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1. Ion Tori, Advection Dominated Accretion Flows 

We will show that the radiative efficiency of an accreting black hole (77 = L/Mc 2 , L is 
the luminosity, M is the mass accretion rate) depends on the magnetic field strength. An 
investigation similar to ours, but with different conclusions, was carried out by Quataert 
(1997). 

It was suggested that at low accretion rates, M <C L Edd /c 2 , a two-temperature plasma 
with Tj ^> T e forms near the Schwarzschild radius of an accreting black hole (Shapiro et 
al 1976, Phinney 1981, Rees et al 1982, Narayan & Yi 1995), and the disk is inflated into 
a torus by the ion pressure. The plasma in the torus is collisionless, and Coulomb energy 
exchange between electrons and ions is small. Viscous heating power must be ~ O.lMc 2 . 
If nearly all the heat is released into ions, radiative efficiencies as low as 77 ~ 10~ 4 are 
obtained, because ions take all the heat into the hole. Electrons radiate efficiently, and 77 
can not be much smaller than 10~ 4 due to adiabatic heating of electrons. Observational 
estimates for L and M exist for a number of objects ( Narayan et al 1997, Fabian & Rees 
1995 ), and some black hole candidates seem to have r] <^ 10~ 4 . 

We show here that if indeed 77 ^ 10~ 4 , magnetic fields in these accretion flows must be 
very weak, (3~ l = B 2 j%i\p <^ 10~ 3 , where p is the gas (ion) pressure. For stronger magnetic 
fields, 77 > 10~ 4 because electrons receive more than 0.1% of the viscous heat, and this part 
of the viscous heat is efficiently radiated. At equipartition fields (f3~ l ~ 1), electrons get 
most of the heat and the radiative efficiency is close to the standard thin disk estimate , 
77 ~ 0.1. Theoretical models of black hole accretion that assume both (3' 1 ~ 1 and 77 <^ 10~ 4 
should be abandoned. 

The reason why stronger magnetic fields lead to an increased fraction of viscous energy 
going into electrons (and correspondingly to larger radiative efficiency) is roughly as follows. 
At strong fields, turbulence in the torus is dominated by Alfven waves which are only 
weakly damped at the length scales ~ p (the ion Larmor radius), where damping by ions 
is maximal. Turbulent energy flux gets through the ion damping barrier at scales ~ p and 
cascades to smaller scales. At length scales <C p, ions drop out of dynamics because of the 
Larmor circle averaging, and the energy flux heats electrons. 

The constraint fi^ 1 <; 10~ 3 (for 77 <^ 10~ 4 ) is of interest for modeling the spectrum of 
an accreting black hole. The radiation spectrum is determined by the synchrotron emission 
and its Comptonization; the synchrotron peak frequency ~ eB/mc changes by a factor of 
30 when our estimate of is changed from 1 (which is the standard assumption) to 10~ 3 .[] 



The peak frequency also depends on the electron temperature. This dependence is non-trivial because 
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2. Electron and Ion Heating 

Unless the plasma somehow leaves precisely all of its angular momentum behind by 
the time it reaches the hole, the ion torus will be rotating at about the Keplerian speed. To 
accrete, the plasma has to transport its angular momentum outwards (in a frame moving 
radially with the plasma), and this means that hydromagnetic turbulence should develop 
in the torus. Whatever the nature of turbulence and the resulting turbulent viscosity are, 
viscous heat will be generated at a rate ~ O.lMc 2 (e.g., Shapiro & Teukolsky 1983). The 
question is into which plasma species is this heat deposited. The ions will take their share 
into the hole, while most of the electrons' share will be radiated. Narayan et al (1997) 
and Mahadevan (1997) suggested that the ratio of electron to ion heating, 5, is given by 
the mass ratio, 5 ~ m e /mi ~ 10 -3 . We will show that such small ratios of electron to ion 
heating are possible only at very weak magnetic fields. 

There are two particle heating mechanisms in a magnetized collisionless plasma - 
Landau heating by electric fields parallel to the local magnetic field and Landau heating by 
the time varying magnetic field (by the curl of the electric field parallel to the local magnetic 
field, e.g. Stix (1992)). In the advection dominated regime (Narayan & Yi 1995) near the 
hole, thermal velocities of electrons and ions are v e ~ Vi ~ c, yet the 7-factors are not large, 
so that the plasma is roughly non-relativistic. At comparable thermal velocities Landau 
heating by parallel electric fields is proportional to Q ~ e 2 /m, and hence it is dominated 
by electrons, the ratio of electron to ion heating being 5 ~ rrii/m e . For magnetic heating, 
the charge e should be replaced by the magnetic moment, Q ~ /i 2 /m, where [i = mv 2 /2B. 
The heating ratio is 5 ~ m e /mj, as Narayan and coworkers usually assume. Electron to 
ion heating ratio will be much greater than if heating by the parallel electric fields is 
important. Ion heating is also suppressed, and correspondingly electron to ion heating ratio 
is increased, at small length scales when averaging of the Lorentz force over the ion Larmor 
circle becomes important. 

To calculate the ratio of electron to ion heating, we must follow the turbulent energy 
as it cascades from the largest energy containing scales to smaller scales and estimate the 
strengths of the two damping mechanisms for all scales. Since random magnetic fields 
limit the mean free path of plasma species to corresponding Larmor radii, hydrodynamic 
description is applicable at large scales, and Kolmogorov energy cascade develops. At 
smaller scales, turbulence is dominated by Alfven waves, because the nonlinear frequency 
(e/c 2 ) 1 / 3 (e is the energy flux in momentum space, units cm 2 /s 3 ) becomes smaller than the 
linear Alfven frequency VAk (va is the Alfven speed). Following Goldreich & Sridhar (1995) 



the synchrotron emission is self-absorbed and comes from a range of distances from the hole. 
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( Appendix A), we assume that for k > e/v\, when Alfven frequency is formally larger than 
the nonlinear frequency, turbulence develops anisotropy, with filamentary perturbations 
aligned with the local magnetic field, k± ^> k\\. Degree of anisotropy is determined by a 
requirement that the spectrum be critically balanced, meaning that hydrodynamic and 
Alfven frequencies are comparable , (e/c^) 1 / 3 ~ VAk\\. 

In the critically balanced Alfven turbulence regime, it takes about t ~ l/i^fey for a 
turbulent eddy to cascade to twice its wavenumber. Let 7^ be the Landau damping rate 
at a given scale (calculated in Appendix B). The fraction of energy flux intercepted at 
wavenumber scales k is of order 7fc/i>A&|| - the rest of energy goes to scales of order ~ 2k. 
The ratio 7/cj (damping to real frequency for linear waves) is all we need. This ratio 
determines at what scales the energy flux is damped, and knowing this it is straightforward 
to predict which plasma species is heated. Figure 1 shows 7/0; as a function of the 
perpendicular wavenumber. As seen from the figure, significant fraction of the energy is 
damped near the ion Larmor radius, k±p ~ 1, when magnetic fields are weak, < 0.1 

For kj_p <^ 1 the main damping mechanism of Alfven waves is magnetic heating of ions. 
The damping rate is given by (Appendix B, for f3~ l <C 1) 

7 = W7^^ p)2 - (1) 

U) o2a/7T Va 

At kj_p ^ 1, Alfven waves gradually modify into whistlers (Appendices B and C), the ions 
still take most of the heat load, and for f3~ l <C 1, 

1 = 1 ■ (2) 

According to the last formula, ion damping becomes small at k±p ~ few. Turbulent energy 
that gets through the ion damping barrier at kj_p ~ 1 is only weakly damped until it reaches 
very small scales to be damped by electrons (see Appendix C for a discussion of whistler 
turbulence). Ions drop out of the picture at k±p ^ few because of Larmor circle averaging. 
The ion heating maximum at k±p ~ 1 determines the ratio of electron to ion heating. 

At equipartition magnetic fields (f3~ l = 1), the maximal value of ion damping is 
7/cj ~ 0.02 at k±p ~ 1 and much smaller at other scales. Most of the viscous heat gets over 
the ion damping barrier and is ultimately damped in electrons. Figure 2 shows the turbulent 
energy flow in the wavenumber space (arrows are drawn by hand!). The corresponding 
radiative efficiency is 77 ~ 0.1, similar to thin accretion disk models. 

The very small radiative efficiency assumed in advection dominated accretion models 
, 7] ~ 10~ 4 , is possible only if a very small fraction, 5 <; 10~ 3 , of the Alfven wave energy 
at k±p < 1 is transmitted to whistlers at kj_p ~ few. Whistlers can propagate at k±p as 
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small as ~ 1. For (3 > 10 3 , Alfven waves are not fully damped at k±p xs 0.5 and will excite 
whistlers by three-wave interactions. Thus, ^ 10~ 3 is necessary for rj <^ 10~ 4 . 

Fortunately, electron parameters (v e and m e ) do not enter Eqs. (1) and (2). Electron 
thermal velocity and 7-factor are determined by radiation and are not reliably predicted by 
the advection dominated (or any other) accretion model. However, the only parameter that 
matters for our calculation is - an input parameter. Our predictions are not limited to 
any particular model of accretion, as long as Tj ^> T e and va < Vi < v e . 

3. Discussion 

A minimal model of turbulence was assumed in this work - no secondary instabilities, 
no shocks, just a simple Kolmogorov cascade with Landau damping. This model predicts 
substantial heating of electrons and large radiative efficiencies if magnetic fields are close to 
equipartition. Small radiative efficiency is possible if nearly all the viscous heat is absorbed 
by ions. For this, magnetic fields must be very weak, j3~ l <^ 1CT 3 . 

Incidentally, we can argue that low radiative efficiency is possible only at very 
weak magnetic fields if we assume that Blandford & Znajek (1977) mechanism is at 
work. Poloidal magnetic field of order B v near a maximally rotating black hole gives 
a luminosity L ~ B 2 R 2 c. We get a radiative efficiency rj ~ {B p / 'B) 2 a~ l , where a 
is the Shakura-Sunayev parameter. Arbitrarily assuming B p /B ~ 0.1, a ~ 0.1, we get 
ijbz ~ 0.L3 -1 , so that (3' 1 < 10" 3 is required for rj BZ < 10" 4 . 

I thank John Bahcall and Martin Rees for suggesting the problem, and Peter Goldreich, 
Pawan Kumar and Insu Yi for valuable discussions. This work was supported by NSF 
PHY-9513835. 



A. Strong Turbulence of Alfven Waves According to Goldreich &i Sridhar 

This Appendix is based on Goldreich & Sridhar (1995). Suppose that e is the energy 
flux in the wavenumber space (energy pumping rate, units cm 2 /s 3 ). Isotropic Kolmogorov 
turbulence is characterized by a nonlinear frequency ~ (ek 2 ) 1 ^ 3 at wavenumbers ~ k. If 
there is a large scale magnetic field with Alfven speed va, another characteristic frequency 
of the turbulence is the linear frequency i>a&||. If the turbulence were isotropic, the 
hydrodynamic nonlinear frequency would be much smaller than the linear Alfven frequency 
at small enough scales. According to Goldreich & Sridhar this does not happen. Instead, 
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turbulent energy propagates to smaller scales within a narrow cone in the /c-space given by 
k\\ & (^ 3 /vA)k'± 3 - In other words, turbulence is dominated by elongated structures, that 
are aligned with the local magnetic field. This could be roughly explained as follows. 

Suppose, that a "turbulent eddy" leaves the cone. On the outside of the cone, 
turbulence is just a collection of weakly nonlinear waves. Waves propagating parallel and 
antiparallel to the local magnetic field follow the field lines perturbed by each other, and 
thus cascade in k±. Turbulent energy from outside of the cone will cascade in k± without 
cascading in k\\ and will thus enter the cone. 

Turbulence fills the entire cone, because within the cone the presence of the large-scale 
magnetic field is not "felt", i^/cy < (e/c 2 ) 1 / 3 . Turbulent eddies will try to isotropize 
themselves, filling as much of the /c-space as possible. 

B. Landau Damping of Alfven Waves and Whistlers 

We calculated Landau damping of Alfven waves (whistlers at kj_p ^ 1) numerically 
using a full dielectric permittivity tensor of a magnetized Maxwellian plasma (e.g. 
Mikhailovskii 1967). Results are shown in Fig.l, where the ratio of imaginary to real 
parts of the eigenfrequency is plotted versus dimensionless perpendicular wavenumber k±p. 
The calculation was carried out for Maxwellian electrons and ions with thermal velocities 
v e /c = 0.66 and Vi/c = 0.33 which are close to the plasma parameters of a two-temperature 
advection dominated accretion model of Narayan & Yi (1995). We assumed k\\/k± = 10~ 2 , 
and checked that does not depend on k\\/k± as long as k\\/k± <C 1. 

The four curves correspond to four values of Damping by ions is maximal 

at k±p ~ 1, where damping by electrons is negligible. For strong fields, ~ 1, the 
damping at k±p ~ 1 is small. The damping increases with decreasing /3' 1 . At ~ 10~ 2 , 
the damping is so strong that Alfven waves do not propagate in a narrow region near 
k±p ~ 0.85. The region of forbidden propagation is narrow for f3~ l ~ 10~ 2 and widens for 

/r 1 ~ 10- 3 . 

Two remarks are in order. Electrons and ions that resonate with Alfven waves and 
whistlers belong to the bulk of the distribution functions. Deviations from Maxwellian 
distribution, which may be expected in our collisionless plasma, will not strongly affect 
the results. Also, linear and nonlinear Landau damping rates are similar, because linear 
frequencies and nonlinear frequency shifts are of the same order for a critically balanced 
turbulence. 
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Most of the numerical results shown in Fig.l can be understood analytically. As 
discussed in §2, analytical expressions (1) and (2) are of importance - they imply that our 
results are model independent. For k± ^> k\\, va <^ v i <; v e <; c, the following approximate 
expressions for the permittivity tensor (Mikhailovskii 1967) give a good accuracy. For a 
wavenumber in the xz-plane, with z along the local magnetic field, 

c 2 I -he- (nr , 

v\ z 

e xy , ig = i ^ 1 L {1 - {Io - Il)e -. } , (B2) 

e X z = 0, (B3) 
e yy = e 2 = i^W{v/v^) 2 ^{h-h)e-\ (B4) 

II ^ A ^ 

e yz = 0, (B5) 

Here v — uj/k\\ is the phase speed, z = (/cj_p) 2 /2, / = I(z) are modified Bessel functions, 
and 

W(x) = e~ x \l + f die 1 ' 2 ). (B7) 

Permittivities ei and e 2 are due to ions, e 3 is due to electrons, and g is due to both species. 
The dispersion relation for linear waves follows from Maxwell equations, 

c 2 

detfe,- + —(kikj - k 2 5ij)) = 0. (B8) 

For kj_p 10, the parallel permittivity e 3 is very large, and both damping and real 
frequency are accurately predicted by a two dimensional approximation 



c 2 kl, c 2 k 2 . 



^-—r)^-—r)-9=^ (B9) 

For k±p <C 1, we have e±e2 3> g 2 , and Eq. (B9) reduces to e± = c 2 /v 2 , which gives the 
Alfven wave dispersion law 

v = v A . (BIO) 

For kj_p > 1, we have c 2 kj/u 2 > ei, c 2 k 2 /u 2 > e 2 , and Eq. (B9) gives c 2 k\\k/uu 2 = g, that 
is whistlers 

v = {v 2 A /v t )pk. (Bll) 
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Corrections to the above dispersion laws (within the 2D approximation, Eq. (B9)) give 
the damping rates for Alfven waves and whistlers. For Alfven waves, in the limit va <C i>j 
when e 2 is pure imaginary, we have approximately 

c 2 k 2 

(ei - —£-)e 2 = g 2 , (B12) 

UJ Z 

giving 

1 = (B13) 
uj 2eie 2 

that is 

- = WKir^pr- (B14) 

uj S2y/ir va 

For whistlers at kj_p ^ 1, we have approximately 

^(^- £2 )= 9 2 , (B15) 

UJ A UJ Z 

giving 

« = 2(i^- <B16) 

At even larger k±p electron damping becomes important. The ratio of electron to ion 
damping can be estimated (for va <C f») as 



\E Z \ 



Am\ zp 12 ' 

e 2 l-C'j/l 



(B17) 



Here the wave electric fields should be obtained from the system of linear equations whose 
determinant is given in Eq. (B8). Two of these equations can be approximately written as 

°-^E x = e 3 E z , (B18) 



c 2 k 2 ± 



E y = igE x . (B19) 



Now Eq. (B17) gives 

6 = ^™{k ± p)\ (B20) 
2 rriiVi 

showing that electron heating dominates for k±p ^10. At these small scales, ions are 
not important either for real part of the frequency or its imaginary part (heating). For 
perturbations with perpendicular length scales <; p/10, ions can be treated as a constant 
positive background for electrons and electromagnetic fields. 
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C. EHD and Whistler Turbulence 

We should make sure that turbulent energy will keep cascading to smaller scales after 
Alfven waves are transformed into whistlers at length scales ~ p. Like collisionless Alfven 
waves, collisionless whistlers can be described hydrodynamically. Hydrodynamic equations 
valid at length scales ^ p can be obtained as follows. Electrons move freely along magnetic 
field lines, and En = 0. In the perpendicular direction, electrons move with the ExB-drift 
velocity, vj_ = cE x B/I? 2 . These two equations give 

E + -vxB = 0. (CI) 

c 

Neglecting displacement currents, 

Att 

V x B = nev. (C2) 

c 

From Eqs. (CI) and (C2), we obtain the Electron Hydrodynamics (EHD) equation 
(Kingsep, Chukbar, & Yan'kov 1990) 

d t B = — ^— V x (B x V x B). (C3) 
Airne 

Linear waves in Eq. (C3) are whistlers, their dispersion law can be presented as Eq. (Bll). 
Turbulent cascade in Eq. (C3) was discussed by Kingsep, Chukbar, & Yan'kov (1990) and 
by Goldreich & Reisenegger (1992). 
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Fig. 1. — Damping rate over real frequency versus k±p for plasma parameters given in 
Appendix B. For k±p < 10, ions are taking most of the heat, with the exception of the 
P -1 = 1 case; the rise of the curve at \ogk±p > 0.5 is due to electron heating. 
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Fig. 2. — Energy flow in the wavenumber space for = 1. Narayan and Yi (1995) 
formulae were used to estimate the parameters. In their notation, black hole mass m ~ 10 6 , 
Shakura-Sunyaev parameter a ~ 0.1, mass accretion rate in Eddington units rh ~ 0.001, 
turbulent energy-containing length scale in units of the Schwarzschild radius L ~ 0.3, large 
scale turbulent velocity ~ 0.1c. 



